Measurement of main parameters of the \J/{2S ) resonance 
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Abstract 
, 

^ A high-precision determination of the main parameters of the t^(2S ) resonance has been performed with the KEDR detector at the 
I VEPP-4M e'^e' collider in three scans of the ifi{2S )-i/'(3770) energy range. Fitting the energy dependence of the multihadron cross 
section in the vicinity of the i//(2S ) we obtained the mass value 



> 



M = 3686.114 ±0.007 + 0.011 'l'^^'^I MeV 
and the product of the electron partial width by the branching fraction into hadrons 

r«, X Si, = 2.233 ± 0.015 ± 0.037 ± 0.020 keV. 



T^lj- The third error quoted is an estimate of the model dependence of the result due to assumptions on the interference effects in the 
cross section of the single-photon e^e" annihilation to hadrons explicitly considered in this work. Implicitly, the same assumptions 

f-->) were employed to obtain the charmonium leptonic width and the absolute branching fractions in many experiments. 

T— I ' Using the result presented and the world average values of the electron and hadron branching fractions, one obtains the electron 
partial width and the total width of the i^{2S ): 

> 



X 



Fee = 2.282 + 0.015 ± 0.038 + 0.021 keV, 
r = 296±2±8 + 3 keV. 



- - These results are consistent with and more than two times more precise than any of the previous experiments. 



1. Introduction 

More than thirty six years passed since the discovery of J/ifi, 
but studies of charmonium states still raise new questions. Re- 
cent progress in charmonium physics requires an improvement 
of the accuracy of the parameters of charmonium states 
This Letter describes a measurement of the ^(2S) meson pa- 
rameters in the KEDR experiment performed during energy 
scans from 3.67 to 3.92 GeV at the VEPP-4M e+e" colUder. 
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For a precision experiment it is essential to state explicitly 
what quantities are measured and how they can be compared 
with results of theoretical studies, therefore we discuss a defini- 
tion of the \I/{2S ) parameters just after a brief description of the 
experiment. The importance of the question has grown since 
the appearance of the work iQ] in which the BES collaboration 
used an original approach to the determination of the j'''- - \ 
resonanceparameters. Its further modification has been used in 
Refs. Sai. 
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2. VEPP-4M colUder and KEDR detector 

VEPP-4M is an e^e" collider j^l designed for high-energy 
physics experiments in the center-of-mass (cm.) energy range 
from 2 to 12 GeV. The peak luminosity in the 2x2 bunches op- 
eration mode is about 2 x 10-"' cm^^s"' in the vicinity of ^{2S ). 
Having a modest luminosity, VEPP-4M is well equipped for 
high-precision measurements of beam energy The instanta- 
neous value of the beam energy can be calibrated using the reso- 
nant depolarization method (RDM) iH^] with the relative accu- 
racy of about 10"^. The results of RDM calibrations must be in- 
terpolated to determine the energy during data taking and the in- 
terpolation accuracy of about 10 keV can be reached lioll . Con- 
tinuous energy monitoringis performed using the infrared light 
Compton backscattering jl ill with the accuracy of the method 
~ 60 keV. 

The KEDR detector ifl?] (Fig.[T]i comprises the vertex de- 
tector (VD), drift chamber (DC), time-of-flight (TOE) system 
of scintillation counters, particle identification system based on 
the aerogel Cherenkov counters, EM calorimeter (liquid kryp- 
ton in the barrel part and Csl crystals in the endcaps), super- 
conducting magnet system and muon system inside the magnet 
yoke. The superconducting solenoid provides a longitudinal 
magnetic field of 0.6 T. The detector is equipped with a scat- 
tered electron tagging system for two-photon studies and some 
applications. The on-line luminosity measurement is provided 
by two independent single bremsstrahlung monitors. The trig- 
ger consists of two hardware levels: the primary trigger (PT) 
and the secondary one (ST) ifisll . The PT operates using sig- 
nals from the TOE counters and fast signals from the Csl and 
LKr calorimeters, whereas the ST uses the normally shaped 
calorimeter signals and the information from VD, DC and TOE 
system. After the readout, a software selection of events is per- 
formed using simplest event characteristics, in particular, the 
number of hits in VD. The upper limit on the number of VD 
tubes hitted is very effective for the machine background sup- 
pression. 

3. Data sample 

In 2004 two scans of the ^{2S )-tfr{3770) energy range were 
carried out with an integrated luminosity of about 0.7 pb"'. In 
2006 the regions of iff{2S) and i^(3770) mesons were scanned 
once again with an integrated luminosity of ^ 1.9 pb The 
combined data sample corresponds to 1 .6 ■ 10^ i//{2S ) produced. 
The data acquisition scenario for tlr{2S ) was similar to that de- 
scribed in Ref. The accuracy of the energy interpolation 
between the RDM calibrations varied from 10 to 30 keV during 
the whole experiment. Our final results on the i/'(3770) param- 
eters are presented in the next letter of this volume. 

4. On definition of J^'^ = 1 resonance parameters 

A resonance with the quantum numbers j'''^ = 1 can be 
treated in some cases not as an unstable particle but just as a 
vacuum polarization phenomenon. Sometimes this causes con- 
fusion in the data analysis as was noted in Ref. fl^. To avoid 




Figure 1: 1 - Vacuum chamber, 2 - Vertex detector, 3 - Drift 
chamber, 4 - Threshold aerogel counters, 5 - ToE-counters, 6 - 
Liquid krypton calorimeter, 7 - Superconducting solenoid, 8 - 
Magnet yoke,9 - Muon tubes, 10 - Csl calorimeter, 11 - Com- 
pensating superconducting coils 



confusion, "bare" and "dressed" or "physical" parameters of a 
resonance must be clearly distinguished. The former do not in- 
clude QED corrections and are used in many theoretical studies, 
the latter include some of them (in particular, the vacuum polar- 
ization) and are published as results of almost all experimental 
papers. 

The physical parameters correspond to the interpretation of 
a 1 resonance as a particle described with a Breit-Wigner am- 
plitude representing its appearance in all orders of the QED per- 
turbations , which is absolutely natural for the strong production 
of a resonance. Eor an electromagnetic production of a reso- 
nance like, e.g., in the process e^e" — > 1^(25 ) — > p'p, this ampli- 
tude interferes with the pure QED amplitude of the e^e^ — » p'p 
transition. As shown below, such an approach allows one to 
avoid the numerical integration in the calculation of radiative 
corrections. 

Let us demonstrate the relation between the bare and physi- 
cal parameters. According to Ref. Iil5il , the cross section of the 
single-photon annihilation can be written in the form 



= J'dx 



crQ{{l-x)s) 

1 -n((i-jc)i)|2 



r(s,x). 



(1) 



where s is the cm. energy squared, Tis, x) is the radiative cor- 
rection function, Yl{s) represents the vacuum polarization oper- 
ator and cro(i) is the Born cross section of the process. One has 



n = IIo H- with the nonresonant IIq 
and the resonant 
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where a is the fine structure constant, Mq, Fq and rfj are the 
"bare" resonance mass, total and electron widths, respectively. 
Eg . (|2]l slightly differs from the expression used in Refs. jig, 
llTll . It corresponds to the simplest resonance cross section pa- 
rameterization 
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o-r(s) = ImnR(5) = 



i27rrl°'ro 



(3) 



For the muon pair production the Born cross section at s » 



4m^ {nift is the muon mass) is just 



4na^ 
^7" 



(4) 



thus the resonance behaviour of the cross section is, in this ap- 
proach, entirely due to the vacuum polarization, which implic- 
itly describes the muonic decay of a resonance. Eqs. ([T]), ^ 
and (HI give the dimuon cross section without separation into 
the continuum, resonant and interference parts. To obtain the 
contribution of the resonance, the continuum part must be sub- 
tracted from the amplitude. It can be done with the identity 
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i-no-n«(,9) 



1 



1 



3r 



,(0) 



1 
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l-Oo (l-no)2 a Mos-M^ + iMF 



Two terms in the right-hand side correspond to the continuum 
amplitude and the resonant one, respectively. The second power 
of the vacuum polarization factor 1/(1 - Do) in the latter can be 
interpreted as the presence of two photons, one at a resonance 
production and the other in its decay. In the resonant amplitude 
both M and f depend on s 
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Mr = MoTo • 



a Mq 1 - Ho 
1 
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Im- 



a Mq 1 - IIo 



In the vicinity of a narrow resonance this dependence is neg- 
ligible, thus the resonant contribution can be described with a 
simple Breit-Wigner amplitude containing the physical param- 
eters M M(Ml) and F * V(MI). 

To obtain the dimuon cross section one has to multiply the 
absolute value squared of the right-hand side of Eq. (|5]l by the 
Born cross section (|4|l. The resonant part of the cross section 
is proportional to F^^* squared which appears instead of F^^^F^^' 
due to the lepton universality in QED. The factor 1/(1 - Oq)^ in 
front of the resonance amplitude converts the square of the bare 
vf] to the square of the physical partial width 



-(0) 



Tee = 



(7) 



recommended to use by Particle Data Group since the work ifisl] 
appeared. 

For electromagnetic decays of a resonance to hadrons, the 
only difference with the dimuon case is the factor R in the cross 



section, where R is the hadron-to-muon cross section ratio off 
the resonance peak. For a strong decay with the partial width 

is) 

Fq the Born cross section is 



{s-Ml)^+MlTl 



(8) 



In this case the identity (|5]) is not required, the direct substi- 
tution of (|2]i and dS) in ([TJ leads to the same definition of the 
physical mass, total width and leptonic width. The equivalent 
definition of the physical mass in the hadronic channel is given 
in Ref. lfl9l] . The physical value of the partial width F*'' is iden- 
tical to the bare one V^^K 

We would like to emphasize that the experimental values of 
a 1 resonance mass, total or leptonic width can not be com- 
pared with the immediate results of potential models or used to 
fit parameters of a potential without either "undressing" of the 
experimental values or "dressing" of the potential model results 
with Eqs. (|6]l and ©. The differences between dressed and bare 
masses are about 1.2 and 0.5 MeV for the Jlifj and iff{2S), re- 
spectively. The corresponding differences for the total widths 
are about 23 and 10 keV (ImFIo ~ -a iR+2)/3 with the R ratio 
about 2.2). 

Unlike the works 12. S, we consistently use the physical 
parameters and treat equally strong and electromagnetic decays 
of the t//{2S). The difference in the approaches is discussed in 
more detail in Section 15.31 It is worth noting that we do not 
suggest any new approach, but just follow the one employed in 
most measurements of heavy quarkonium parameters though, 
as far as we know, its relation to Eq. ([TJ in the hadronic and 
leptonic channels was not rigorously considered until recently. 

5. Cross section calculation 

The cross section formulae given below in this section con- 
tain the 4i{2S) total width F and the electronic width F^e, the 
determination of which are the goal of our analysis. We fix 
these parameters in the cross section fit, but use the iteration 
procedure to obtain the final results. The values of the i/'(25) 
electron and hadron branching fractions, required to recalculate 
Tee X 'Bit to Vee and F, are fixed at the world averages. The sys- 
tematic uncertainties due to such an approach are discussed in 
SectionlT4l 



5.1. Multihadron cross section 

Below we present the results of the paper [Eo'l published 
soon after Jj^ discovery in the updated interpretation. Some 
details of the analytical calculations and numerical checks can 
be found in Ref. 1I21I1 . 

Using the physical values of the parameters, for strong de- 
cays of \I/{2S ) one reduces Eq. ([1]) to 



\2nTeeT^ 



{W-{\-x)- 



■ M^y + M2F2 
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where W - ^/s is the total colUsion energy, T, F^e and fJ*' are 
the total and partial widths of the ip{2S) meson, and M is its 
mass. 

Taking into account the resonance-continuum interference 
and performing the integration over x with a simplified version 
of T^ix, s) one obtains 



TM 



ImfiW) 



A Re — - 

3W l-Hn 



2TM 
TM 

In 



1 

1 H r arctan - . . _ 



(10) 
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where Ho is the vacuum polarization operator with the i/f(25 ) 
contribution excluded. The f/, parameter includes both strong 
and electromagnetic decays and some contribution of interfer- 
ence effects which is discussed in the next subsection. 

The first square bracket in Eq. ( fTOt corresponds to radiation 
of soft photons, while the second one represents hard photon 
corrections. The A parameter introduced in Ref. 1I20I1 charac- 
terizes the strength of the interference effect in the multihadron 
cross section and equals 1 for the dimuon cross section. 

The correction 5sf follows from the structure function ap- 
proach of Ref. I15II : 



3 a In- 



2 1 
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2/37 1 
^^'96 - 12 - 36^" 



(11) 



40-/, W 1\ 

In— -T , (12) 
71 \ nie 2] 

nie is the electron mass and the function / is defined with 



/(W) = 



sinnp \m2 - W2 - iMT 



(13) 



The presentation of the soft photon integrals in form of the real 
and imaginary parts of the function / is more transparent than 
that of Ref. 

Despite a simplification of 'Fix, s), not too far from the 
i//(2S ) peak the resonant part of Eq. ( fTOt reproduces the results 
obtained by the numerical integration of the complete formula 
with an accuracy better than 0.1%. 

The resonant part of Eq. ( fTOb is proportional to the Ffgf /,/r 
combination which is a product of the partial width and the 
branching fraction: Fee x Si, - F/, x See- Since our final re- 
sult is Fee X S/,, let us consider the relation of f with the true 
Ff,, which is a sum of hadronic partial widths. 

5.2. Interference effects in total multihadron cross section 

Considering charmonium decays at the parton level, one 
deals with the gluonic if/ ggig/y) and electromagnetic ifr 



y* — > qq modes. Treating quarks and gluons as final decay prod- 
ucts, one obtains that gluonic modes do not interfere with the 
continuum e^e^ -^y* ^ qq process while those of electromag- 
netic origin do with the interference phase equal to that of the 
dimuon decay. In this case F/, does not differ from a sum of the 
hadronic partial widths F/, - Fgg^gy) + Fqq with Fqq - R F^g and 



RSe. 



(14) 



where See and S/, denote the electron and hadron branching 
fractions, respectively. The real situation is much more compli- 
cated. 

For an exclusive hadronic mode m at a given point in the 
decay product phase space the ampUtude e^e" — > m can be 
written as 



A J®) 



12n la am(Q) 
[3 l-Ho 



a„m M yfRjfg + al:X&) e'-^"' M yfrj^: 
M2 - W2 - IMF 



(15) 



where R,„ is the mode contribution to R, F,„ represents the 
contribution of the strong interaction to the partial width and 
<p,n is its phase relative to the electromagnetic contribution F)^ = 
Rm Fee, the real functions a are normalized with J a^,(®) d@ = 
1 . In general, the phase 0,„ depends on 0. The numerator of the 
last term of Eq.lfTsTi is proportional to the decay amplitude, the 
partial width is 



RmTee + F^f + 2 ^R„,FeeF)^' <COS 0,„)0 
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where the angle brackets denote averaging over the product 
phase space: { x(0) )0 = / a(0)a*''x(0) d®. To obtain the ex- 
clusive mode cross section, the following replacement must be 
done in the expression ( fTOt : 



F,n — > F/, 
I 

Re 



1-no 



i-no ' 



where the latter replacement follows from comparison of the 
interference term corresponding to Eq. ( fTSb and that of Eq. (fTOl l. 
Performing them and summing over all hadronic modes one 
obtains the expressions for f /, and A:. 



■Fh X 



1 -H 



2a 



3(l-Reno)S,, VS, 
(here F/, = 2 F,„, Imllo is neglected) 



~R I 



(18) 



A = 




Z7„,S^" (cos<^„)0 , (19) 



where b,„ = R„ /R is a branching fraction of the continuum pro- 
cess and — r*^Vr. Below, the sums containing {smcf>„,)Q 
and {cos(P,„)q are referred to as Ssin and Ecos, respectively. The 
parton level results are reproduced by Ea.dTSll and ( fT9] l if Ssi„ 
and Scos can be neglected. For a hypothetical heavy charmo- 
nium decaying to light hadrons, both the values of (cos (pm}e 
and (sin (f>m)@ averaged over the product phase space tend to 
zero due to the different configuration of jets in electromag- 
netic and strong decays. For the real J/i// and 4r{2S) one has to 
rely on the absence of the phase correlations in different decays. 
For the quasi-two-body decays such correlations are expected 
(Ref. 12311 and references therein) but their branching fractions 
are small ll24ll . 

If the sum S,,,, is not negligible, the method of the resonant 
cross section determination employed in this and many other 
experiments becomes inaccurate and ambiguous because of the 
well-known ambiguity in the partial width determination which 
takes place for each individual mode. Indeed, a fit of the mode 
cross section cr„, gives the values of f ,„ and cos 0„, but leaves 
unknown a sign of sin (p^ required to obtain F,,,. 

Equating all (sin 4>i„)@ in Eq. ( fTST i to unity, one sets an upper 
Umit on the inaccuracy of the hadronic partial width AF;, and 
the resonant cross section at the peak Acr^'^'*(Wpeak) used for the 
determination of the branching fractions: 



peak. 



Cr^^^(Wpeak) 



Th 3Sh V See Zj V 



(20) 



For tfr{2S ) the sum in the right part < 1 - S j;^+x, thus Eq. ( l20l l 
gives about 4%. A better estimate employing the A value ob- 
tained with the cross section fit is discussed in Section |8] Until 
this section we omit the tilde mark wherever possible, thus F/, 
and Si, should be read as f /, and S/,. 

Eqs. (fTSl l and ( fT9] l show that the correct account of inter- 
ference effects is essential for a determination of the i/r(3770) 
parameters due to its small value of See and large branching 
fraction to D mesons, nevertheless, it was ignored in most of 
published analyses. The interference effect is crucial for a de- 
termination of the non-DD decay fraction of i/'(3770) as was 
emphasized in Ref. lEsIl . 

5.3. Observed multihadmn cross section 

The multihadron cross section observed experimentally in 
the vicinity of i//(2S ) can be parameterized as follows: 



+ £rr<nt(^) + t^co„t(W). 



(21) 



Here e^(2S} and e^r are the detection efficiencies and their de- 
pendence on W can be neglected. The continuum t^t" cross 
section cr^^ is included according to Ref. ll26ll to extend the va- 



Udity of dlB beyond the i//{2S )-iA(3770) region. 

For the iJ/(2S ) cross section (fTOl i includes the t contribution 
and the A parameter is modified properly: 



RSee 



Srr R^Srr 

e^as) \ R Sh 



(22) 



with = (Tt+t- /(Tfi+ft- ~ 0.39. The reduction of t^t^ detec- 
tion efficiency of about 0.3 compared to the multihadron one is 
accounted explicitly. 

In Eq. (ISTT i this cross section is folded with the distribution 
over the total collision energy which is assumed to be quasi- 
Gaussian with an energy spread crw- 



G(W,W ) - — = exp 

\ 2cr^ 



y/27rcrw \ 
The preexponential factor can be written as 

I + aA + b 



l+ba% 



(23) 



(24) 



It is due to various accelerator effects such as the j6-function 
chromaticity. We fix a = = in our fit and consider the cor- 
responding systematic errors in the Tee x Sh product in Sec- 
tion [73] The presence of this factor and other accelerator- and 
detector-related effects yield too large systematic uncertainties 
when the total width parameter F is left floating in the fit at 
o-iv > 5 X F/2. 

Since the interference effect is included in c^^^^e^, the con- 
tinuum cross section is a smooth function, which with the suf- 
ficient accuracy can be parameterized with 



O-conAW) = Cro 



Wo 

w 



(25) 



where ctq is the value of the background cross section at a fixed 
energy Wo below the tl/(2S) peak. 

In contrast with the commonly used interpretation of the 
cross section as a sum of the resonant, continuum and interfer- 
ence parts employed, in particular, in Ref. [2], in the works 

it is interpreted as a sum of the two parts only: the cross 
section of the resonance and the "nonresonant" one. The lat- 
ter is calculated using the full vacuum polarization operator 
Uf)+n^(2S) (Eq. (10) of Ref. f^). The two approaches are equiv- 
alent provided that the bare parameters enter the Tl^(2S) and the 
electromagnetic contribution is excluded from the cross section 
of the resonance (Eq. (3) of Ref. [3]). The full vacuum po- 
larization operator describes not only the interference, but the 
electromagnetic decays as well. If it is not done, the electron 
width extracted from the cross section fit would have a negative 
bias of about R-See/Sh ~ 0.018. 

5.4. Observed e^e^ cross section 

Bhabha scattering events detected in the calorimeter were 
employed for luminosity measurements (see Section 16.31 for 
more detail). For the large angle Bhabha scattering the contri- 
bution of if'{2S) decays is not negligible. The differential e^e^ 
cross section can be calculated with 



[dnj "1^/q)qed ^ 

1 ( 9 F^ 
--{ -^(1 +cos2 0) 
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The first term represents the QED cross section calculated with 



dIV 
TV 



the Monte Carlo technique 1281 12911 . The second (resonance) 
and the third (interference) terms have been obtained in ll20ll . 
The corrections to the latter are not calculated precisely, but that 
does not limit the accuracy of the published results. Parameters 
j6 and / are defined by Eqs. (fT2l i and ( fTSl l. respectively. 



6. Data analysis 

6.1. Monte Carlo simulation 

The simulation of the experiment was performed in the frame 
of the GEAA^r package, version 3.21 

The if/{2S) decays and the continuum multihadron events 
were generated with the tuned version of the BES generator 13 lil 
based on the JETSET 7.4 code 13211 . At the original parameter 
settings of the BES generator the difference of the simulated 
charged multiplicity and that observed experimentally exceeds 
1%, thus the bias in the detection efficiency up to 2% is ex- 
pected. The procedure of the parameter tuning is discussed in 
detail below in Section 17.21 The decay tables were updated 
according to the recent PDG edition ll24ll . The results are pre- 
sented in Fig.|2] where the most important event characteristics 
obtained in the experiment are compared with those in simula- 
tion. Good agreement is observed. 

The detection efficiency for r^r" events was obtained using 
the KORALB event generator ll33ll . Bhabha events required for 
the precise luminosity determination were simulated using the 
BHWIDE generator ii. 



6.2. Trigger efficiency and event selection 

To reduce systematic errors due to trigger instabilities and 
uncertainties in the hardware thresholds, both experimental and 
simulated events pass through the software event filter during 
the offline analysis. It recalculates the FT and ST decisions with 
tighter conditions using a digitized response of the detector sub- 
systems. To suppress the machine background to an acceptable 
level, the following PT conditions were used by OR: 

• signals from > 2 barrel scintillation counters , 

• signal from the LKr calorimeter in the scan of 2006 , 

• coinciding signals of two Csl endcaps. 

Signals from two particles with the angular separation > 20 
degrees satisfy the ST conditions which are rather complicated. 
The MC simulation yields the trigger efficiency of about 0.96 
for t//{2S) decays. Because of a problem with electronics, the 
LKr calorimeter was not used in the analysis of 2004 data and 
that decreased the trigger efficiency to 0.91. 

The performance of the detector subsystems and the ma- 
chine background conditions were very different in 2004 and 
2006, so that the selection criteria are also different. 

2004 data, first and second scans: 

• > 3 charged tracks , 

• > 2 charged tracks from a common vertex in the interac- 
tion region ( p < 7 mm, |z| < 1 30 mm) , 
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Figure 2: Properties of hadronic events produced in il/(2S) decay. 
Here is the number of events and pj is the transverse mo- 
mentum of a charged track. All distributions are normalized to 
unity. 



• event sphericity S^h > 0.05 . 

Here p and z are the track impact parameters relative to the 
beam axis and z-coordinate of the closest approach point. 
The sphericity parameter is defined as 



3 llfli 
- mm ^ 

2 Zp^ 



(27) 



where summation is performed over all particles of the event 
and the minimum is taken over directions of the axis relative 
to which the transverse momenta prj are calculated. S^h is 
calculated using charged tracks only. The cut on Sch is effi- 
cient for suppression of the e^e^ — > e^e y background, that 
of cosmic rays and some kinds of the machine background, 
though it also suppresses the leptonic modes of the cascade 
decay tf>{2S) — > J I if/ + neutrals (see the low sphericity peak in 
Fig. 13. 
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Figure 3: The muhihadron cross section as a function of the cm. energy for three scans. The curves are the results of the fits. All 
data are corrected for the efficiency, the peak cross section depends on the energy spread cru^. 



2006 data, third scan: 

• > 3 charged tracks or two tracks with the acollinearity 
> 35 degrees, 

• > 2 charged tracks from a common vertex in the interac- 
tion region ( p < 7 mm, |z| < 1 30 mm) , 

• > 1 photons with energy > 100 MeV in the calorimeter, 

• event sphericity Sch > 0.05 . 

• the energy deposited in the calorimeter > 450 MeV . 

Analyzing the third scan we also used the alternative se- 
lection criteria without a tight cut on the sphericity, but with 
additional requirements on the calorimeter response. It allows 
us to check the systematic error due to the sphericity cut. 

For additional suppression of the background induced by 
cosmic rays a veto from the muon system was required in the 
cases when more than two tracks did not cross the interaction 
region or the event arrival time determined by TOF relative to 
the bunch crossing was less than -5 ns or larger than 10 ns. This 
condition was common for all three scans. 

The conditions described above reduce the physical back- 
ground contributions which do not scale with energy like l/sto 
a negligible level, except the tau pair production which we took 
into account explicitly in the observed cross section. The con- 
tribution of the beam-wall and beam-gas events, cosmic events 
and their coincidences was evaluated using data collected with 
the separated beams (about 10% of the full data sample). It 
was about 2% of the observed continuum cross section for the 



third scan and about 0.4% for the first two. The analysis of the 
event vertex distribution along the beam axis confirmed these 
estimates. We did not perform the background subtraction in 
each data point as was done in the work |27] for precise R mea- 
surements. This is not required for the resonance parameter 
determination, the corresponding uncertainties are discussed in 
Section 1731 Simulation of i^(25) decays yields the detection 
efficiencies of 0.63 and 0.72 for the two sets of selection crite- 
ria, respectively. To ensure the detection efficiency stability, all 
electronic channels malfunctioning in some runs during a scan 
were excluded from the analysis of all its runs. 

6.3. Luminosity determination 

The stability and absolute calibration accuracy of the 
bremsstrahlung monitors used for on-line luminosity measure- 
ments (Section 12]) are not sufficient for the precision cross sec- 
tion analysis, thus events of Bhabha scattering were employed 
for the off-line luminosity determination. In 2004 it was pro- 
vided by the endcap Csl calorimeter (the fiducial region 20° < 
e < 32° and 148° < 6* < 160°). In the analysis of the 2006 data 
the LKr calorimeter was employed (40° <9 < 140°) while the 
Csl one served for cross-checks only. 

The criteria for e^e" event selection using the calorimeter 
data are listed below: 

• Two clusters with the energy above 0.25 of the beam en- 
ergy and the angle between them exceeding 165 degrees, 

• The total energy of these two clusters exceeds 1.05 times 
the beam energy. 
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• The calorimeter energy not associated with these two clus- 
ters does not exceed 10% of the total. 

The loose energy cuts were chosen to reduce the influence of the 
calorimeter channels excluded from the analysis as was men- 
tioned above. The tracking system was used only to reject the 
background (e^e^ — > yy and e'^e~ hadrons). The number of 
extra photons was required to be less than two for the additional 
suppression of the latter 

6.4. Fitting procedure 

The collision energy W was assigned to each data acquisi- 
tion run using the interpolated results of the beam energy mea- 
surements and assuming W = 2Ebeam- The runs with close W 
values were joined into points with the luminosity-weighted 
values Wi (i is the point number). 

The numbers of hadronic events A^, and events of Bhabha 
scattering «, observed at the i-th energy point were fitted as a 
function of W using the maximum likelihood method with a 
likelihood function 



obs 



+n°^' In 



+ n. 



(28) 



where a^^p*"*"*' and n^P*"*"*' are the expected (observed) num- 
bers of the hadronic and Bhabha events, respectively. The ex- 
pected numbers of the hadronic and Bhabha scattering events 
were determined as follows: 



here cri^adr and cr^+g- are defined by 



(29) 



CThadr = <2S)(^') + Se^e-(l,adr) ^^^(H^,) , 
CTe+e- - Cr°X-{Wi) + Shadr(e+e-) 0"^(2S)(^;)' 



.obs 



(30) 



where cr^^95,(VV,) and cr°^'^^_{Wi) were calculated according to 
(Tli and by integration of (|26] |. respectively. The detection 
efficiencies entering the formulae were determined separately 
at each point using the run-dependent Monte Carlo simulation. 
The values of the cross-feed selection efficiency Ee*e-(hadr) (the 
probability to select the e^e" — > e^e" event as the hadronic one) 
obtained by MC are about 0.006%, 0.006% and 0.37% for three 
scans, respectively. The corresponding values of Eiwdr(e*e-) are 
0.03%, 0.03% and 0.25%. 

The integrated luminosity L, at the energy point / can be 
derived from the condition (9/(9L, In X = , giving 



Li 



Cr hadr + O- e 



(31) 



Using the likelihood function that takes into account both A^, 
and «, ensures a correct estimation of the statistical uncertainty 
in the fit results. 



Table 1: The main results of the scan fits (statistical errors only 
are presented). 





M, MeV 


r,, X Sh , keV 


P(J(% % 


Scan 1 
Scan 2 
Scan 3 


3686.102 ±0.018 
3686.130 ±0.013 
3686.108 ±0.010 


2.258 ±0.033 
2.229 ± 0.024 
2.226 ± 0.022 


15.8 
29.5 
79.5 



The data of each scan were fitted separately, the free param- 
eters were the ip{2S) mass M, x S/,, the energy spread crw 
and the continuum cross section magnitude ctq. The A parame- 
ter was fixed at the value of 0. 13 according to Eq. (l22l l. The data 
points (cr,, Wi) and the fitted curves are shown in Fig. [3] The 
results of the fits are presented in Table [T] 

The statistical accuracy of the mass values is significantly 
better than that for the three scans performed at VEPP-4M in 
2002 1 10]. The combined analysis is required to take into ac- 
count properly the systematic errors of the six mass values. 
It will be described in a separate paper discussing numerous 
accelerator-related effects. For this reason we omit such a dis- 
cussion below and just present the result and error estimates. 
The analysis of the three values obtained in this work gives 

M = 3686.1 14 ± 0.007 ± 0.01 1 MeV. 

The model dependence of the mass value was estimated float- 
ing the A parameter in the fit. 

7. Discussion of systematic uncertainties in Tee x S/, 

The dominating sources of the systematic uncertainty in the 
Fee X Sh = Fee X f/j/r valuc are discussed in the following 
subsections. The issue of the difference between f/, and the 
sum of hadronic partial widths because of possible correlations 
of interference phases is addressed in the next section. 

7.7. Systematic error of absolute luminosity determination 

The major contributions to the uncertainty of the absolute 
luminosity determination are presented in Table |2l 

Table 2: Systematic uncertainties of the luminosity determination 
in % for three scans. The correlated parts of the uncertainties 
are also presented. The uncertainties for the first and second 
scans are assumed to be fully correlated. 



Source 


Scan 1 


Scan 2 


Scan 3 


Common 


Calorimeter calibration 


0.3 


0.3 


0.7 


0.2 


Calorimeter alignment 


0.3 


0.3 


0.1 


0.1 


Polar angle resolution 


0.8 


0.8 


0.2 




Cross section calculation 


0.5 


0.5 


0.5 


0.4 


Background 


0.1 


0.1 


0.1 


0.1 


MC statistics 


0.1 


0.1 


0.1 




Variation of cuts 


<1.3 


<1.4 


<0.9 


<0.2 


Sut7i in quadrature 


s;1.6 


s;1.7 


s;1.2 


-0.5 
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The uncertainty due to the calorimeter energy calibration 
was estimated by variation of simulation parameters concern- 
ing the GEANT performance, the light yield uniformity in Csl 
crystals, a sensitivity to the energy loss fluctuations between 
LKr calorimeter electrodes etc. Additionally, the Csl individual 
channel calibrations using cosmic rays and Bhabha events were 
compared. 

The LKr calorimeter was aligned with respect to the drift 
chamber using cosmic tracks reconstructed in the DC. Then 
eight modules of the Csl calorimeters were aligned relative to 
the LKr calorimeter using cosmic rays reconstructed in the LKr 
strip system. The beam line position and direction were de- 
termined using the primary vertex distribution of multihadron 
events. The inaccuracy of the alignment resulted in the error of 
the luminosity measurements of about 0.3% for Csl and about 
0.1% for LKr. 

The difference in the polar angle resolutions observed ex- 
perimentally and predicted by MC causes an uncertainty in the 
luminosity measurement, since events migrate into or out of the 
fiducial volume. These errors are 0.8% and 0.2% for the Csl 
and LKj calorimeters, respectively. 

The uncertainty of the theoretical Bhabha cross section was 
estimated comparing the results obtained with the BHWIDE lEsll 
and MCGPJ 12911 event generators. It agrees with the errors 
quoted by the authors. The value of 0.5% quoted in Table |2] 
includes also the accuracy of Eq. (|26] |. 

The background to the Bhabha process from the if){2S ) de- 
cays and reactions e^e^ — > fifi{y) and e^e^ — > jy was estimated 
using MC simulation. It contributes less than 0.3% to the lumi- 
nosity. The resonant part of the background contribution was 
taken into account in the fit (Section [6.41 ). The residual lumi- 
nosity uncertainty due to background does not exceed 0.1%. 

In order to estimate the effect of other possible sources of 
uncertainty, the variation of the cuts was performed within the 
fiducial region in which good agreement between the MC sim- 
ulation and experiment is observed. The cut on the deposited 
energy was varied in the range of 55 - 75% of the cm. en- 
ergy. The cuts on the polar angle were varied in a range much 
larger than the angular resolution, the variation in the Bhabha 
event count reaches 40%. The variations discussed above cor- 
respond to a systematic uncertainty shown in Table |2l These 
eff'ects can originate from the already considered sources and 
statistical fluctuations, nevertheless we included them in the to- 
tal uncertainty to obtain conservative error estimates. 

Finally, we compared an integrated luminosity obtained us- 
ing the LKr and Csl calorimeters in the scan of 2006. The dif- 
ference of about L 1 ± 1 .0% was found which is consistent with 
the estimates in Table |2] 

7.2. Uncertainty due to imperfect simulation of i//(2S) decays 

The imperfect simulation of i/'(25) decays contributes sig- 
nificantly to the Fee X Sh systematic error caused by the detec- 
tion efficiency uncertainty. Eventually, this uncertainty is de- 
termined by the experimental statistics available for the event 
generator tuning and the ability of the latter to reproduce dis- 
tributions of parameters essential for the event selection and 



their correlations. In our case such parameters are the charged 
multiplicity and event sphericity calculated using momenta of 
charged tracks. 

The selection criteria described in Section l672l reject low- 
multiplicity events. The corresponding branching fractions are 
either negligible like for t//{2S) — > tt^tt^ or well measured like 
for if/{2S ) — > J/iJ/ Ji-^ fi^fTyy . Such decays are simulated in the 
BES generator 13111 explicitly. The variation of decay table pa- 
rameters within their errors (Ref. indicates the detection 
efficiency uncertainty of less than or about 0.3%. The dominat- 
ing uncertainty comes from decays of higher multiplicity which 
are simulated using the parton approach and the fragmentation 
model incorporated in the JETSET 7.4 code. 

About 60% of the t//{2S) decays include J/il/ which is at 
the edge or even beyond the region where the JETSET results 
can be trusted, nevertheless it has enough options and parame- 
ters to achieve good agreement with experiment for major event 
characteristics (Fig.|2]i and to estimate the detection eflSciency 
uncertainty. To do this we iterated as follows: 

1 . select a critical option or parameter and modify it using 
an educated guess; 

2. select a complementary parameter and modify it to find 
the value at which the observed charged multiplicity agrees 
with experiment; 

3. calculate the detection efficiency and compare it with pre- 
vious results to estimate the uncertainty. 

In addition to the charged multiplicity, the observed distribu- 
tions in the sphericity parameter, the invariant mass of the pairs 
of the opposite signs and the inclusive momentum spectrum 
were controlled. The versions of the simulation obviously con- 
tradicting to experiment were rejected. The results are pre- 
sented in Tableland illustrated by Fig.|4l 
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Figure 4: Detection efficiency vs. charged multiplicity for differ- 
ent versions of the \^{2S') decay simulation. The fit lines cor- 
respond to variation of one of the selected parameters. The 
shadow box corresponds to the statistical error of the charged 
multiplicity. The statistical error of the efficiency ~ 0.001 is not 
shown. 

The following JETSET options were studied: 
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Table 3: Comparison of different versions of the MC simulation of ifr(2S) decays. 



Version 


JETSET modifications 


Charged muhipHcity 


Efficiency, % 


LUND fragmentation function 


Probability of vector meson formation 




Pw 


a-pj. , GeV 






1 


0.50 


0.55 


4.1391 ±0.0031 


71.392 ±0.090 


2 


0.50 


0.65 


4.1555 ±0.0031 


72.232 ± 0.090 


Wstop, Wniin, o-pj^ varied 




Wstop , GeV 


WniinGeV 


(Tpj , GeV 






6 


0.47 


0.8 


0.65 


4.1429 ±0.0031 


71.783 ±0.090 


7 


0.52 


0.8 


0.65 


4.1488 ±0.0031 


72.049 ± 0.090 


8 


0.56 


0.8 


0.65 


4.1512 ±0.0031 


72.170 ±0.090 


SW^top, CTpT varied 






o-pr , GeV 






9 


0.17 


0.7 


4.1407 ±0.0031 


71.890 ±0.090 


10 


0.17 


0.65 


4.1552 ±0.0031 


72.430 ± 0.090 


Wmin, o-p^ varied 




W„,i„ , GeV 


o-pr , GeV 






11 


0.8 


0.675 


4.1401 ±0.0031 


72.033 ± 0.090 


12 


0.8 


0.65 


4.1529 ±0.0031 


72.378 ± 0.090 


Switched off parton shower 




o-p.^ , GeV 






13 


0.65 


4.1409 ±0.0031 


72.118 ±0.090 


14 


0.55 


4.1554 ±0.0031 


72.709 ± 0.090 


Field-Feynman fragmentation function 




W„op , GeV 


o-pr , GeV 






3 


0.62 


0.58 


4.1372 ±0.0031 


71.475 ±0.090 


4 


0.62 


0.50 


4.1491 ±0.0031 


71.981 ±0.090 


5 


0.62 


0.43 


4.1650 ±0.0031 


72.755 ± 0.090 



1 . LUND fragmentation function, parton showers are on; 

2. LUND fragmentation function, parton showers are off; 

3. Field-Feynman fragmentation function, parton showers 
are on; 

4. independent fragmentation, all momentum conservation 
options. 

No acceptable versions were obtained in the latter case. For 
the LUND fragmentation function the influence of a few most 
critical parameters was investigated. 

The main parameters controlling the fragmentation process 
in JETSET are PARJ(32), PARJ(33) and PARJ(37) (we refer to 
them as Wmin, W^top and ^Wstop, respectively, in Table[3]and be- 
low). The fragmentation of a colour-singlet system formed by 
initial partons or during the parton showering proceeds while 
the energy is greater than Wmin plus quark masses, otherwise a 
pair of hadrons is to be produced. The Wstop parameter serves 
to terminate the fragmentation and produce a final hadron pair 
earlier, taking into account the mass of the last quark pair pro- 
duced. To avoid artifacts in the hadron momentum spectrum, 
the stopping point energy is smeared using the SW^top parame- 
ter (20% by default). The transverse momentum of quarks ap- 
pearing during the fragmentation is controlled by the parameter 



PARJ(21) (cTp,). The parameter PARJ(1 1) (Py) determines the 
probability that a light meson formed in the fragmentation has 
spin 1. The default values of the parameters W,nin, Wstop, o-pr 
and Pv in JETSET are 1., 0.8, 0.36 GeV and 0.5, respectively. 
The corresponding values set in the BES generator are 1., 0.6, 
0.5 GeV and 0.6. 

Figure |4] shows the values of the detection efficiency and 
observed charged multiplicity for the sets of options and pa- 
rameters listed in Table [3] The dashed line corresponds to the 
charged multiplicity observed in experiment with the selection 
criteria of the 2006 scan, (A^exp.) = 4.1494 ± 0.0054. The max- 
imum difference between the detection efficiency values at the 
point of the multiplicity agreement is 0.55% with the central 
value of about 72.19%. This is presented with the segment be- 
tween the arrows and its middle point in Fig.H] The experimen- 
tal and simulated multiplicities have statistical errors, therefore 
the segment of the agreement transforms to the rectangle shown 
with the shadow box. Taking its angles we found the confidence 
interval of (72.16 ± 0.57)% and from that derived that the rela- 
tive uncertainty of the detection efficiency due to ambiguity in 
the choice of the JETSET parameter set is about 0.8%. A very 
similar central value of 72.1 1% was obtained by averaging the 
fourteen efficiencies from Table [3] with the weights inversely 
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proportional to the sum of for the four event characteristics 
under control. 

There is a systematic error in the observed multiplicity re- 
lated to the track reconstruction efficiency, which is not exactly 
the same for the experimental data and simulation. The dif- 
ference was studied using Bhabha events and low-momentum 
cosmic tracks and the appropriate correction was introduced in 
the detector simulation. However, the inaccuracy of the correc- 
tion increases the shadow box size in Fig.|4]thus increasing the 
detection efficiency uncertainty. For the first two scans the ef- 
fect is about 0.4% and it grows up to 0.7% for the third scan 
because of some problems in the drift chamber 

We repeated the procedure described above with the alter- 
native set of event selection criteria and obtained a slightly dif- 
ferent uncertainty estimate. To account for this we introduced 
an additional error of 0.3%. 

The contributions to the detection efficiency uncertainty due 
to imperfect simulation of i^(25' ) decays are summarized in Ta- 
ble S 

73. Detector- and accelerator-related uncertainties in Y^e xSh 
The major sources of the systematic uncertainty in the Fee x 
!Bh product are listed in Table |5] 

The systematic error related to the efficiency of the track re- 
construction was considered in the previous section. The trig- 
ger efficiency uncertainty is mainly due to the uncertainty of the 
calorimeter thresholds in the secondary trigger. The estimate of 
about 0.2% was obtained varying the threshold in the software 
event filter. The excess of the filter threshold over the hardware 
one in units of the hardware threshold width was varied from 5 
to 6. 

The trigger efficiency for all scans and the event selection 
efficiency for the third one depend on the calorimeter response 
to hadrons. Their nuclear interaction was simulated employing 
the GEISHA iQ and FLUKA Q packages as they are im- 
plemented in GEANT 3.21 Q. Both of them satisfactorily 
reproduce the pion signal in Csl, but in liquid krypton the per- 
formance of GEISHA is much better, thus we determined the 
efficiencies using it and estimated the systematic errors com- 
paring the results obtained with two packages. 

The crosstalks in the vertex detector electronics introduced 
a variation of the detection efficiency up to 1.5% because of the 
cut on the total number of VD hits in the on-line event selection 
software (Section[2]). A code simulating the crosstalks was de- 
veloped and tuned using either events of Bhabha scattering or 
cosmic rays. The residual uncertainty thus determined is about 
0.1h-0.2% depending on the VD voltage. 

The effect of other possible sources of the detector-related 
uncertainty was evaluated by varying the event selection cuts. 
The conditions on the number of tracks were tightened one by 
one for all three scans. For the last one a cut on the energy de- 
posited in the calorimeter was also increased by the most proba- 
ble photon energy. The detection efficiency varied from 0.53 to 
0.63 in the 2004 scans and from 0.70 to 0.79 in the 2006 scan. 
The maximum variations of the Fgg x S/, result are presented in 
Table 13 



The systematic error in the T^e x Si, value due to the beam 
energy determination (Ref. 101) and the data point formation 
(Section|631i was studied for each scan. The relative error does 
not exceed 0.2% except the third scan which includes one point 
with a significant accelerator instability that increases the error 
up to 0.6%. 

The non-Gaussian effects in the total collision energy dis- 
tribution contribute about 0.2% to the r^e x Si, uncertainty. 
Changing the zero value of the a parameter entering the preex- 
ponential factor (l24l i to the value measured with the specific ac- 
celerator technique leads to the mass shift of a few keV ifioll and 
causes a negligible bias in the x Si, value which is related 
to the area under the resonance excitation curve. The quoted 
estimate was obtained by releasing the b parameter 

The bias in the resonance parameters due to the admixture 
of the machine and cosmic background to selected multihadron 
events was evaluated with the controllable increase of this ad- 
mixture. For each data point the background event sample was 
formed containing the events which passed some loose selec- 
tion criterea but were rejected by the multihadron ones. The 
number of multihadrons events was increased by such frac- 
tion of the background sample size, that the data fit yieled 2% 
(0.4%) growth of the continuum cross section for the third (first 
and second) scan in accordance with the backround estimates 
quoted in Section 16.21 The variation of the il/{2S ) mass value 
did not exceed 1 keV, the change of the F^f x Si, product was 
less than 0.1%. 

7.4. Other uncertainties 

In this subsection we discuss the uncertainties related to 
the iterations used to obtain the total and electron width val- 
ues (Section|5]l, the fixation of the interference parameter A en- 
tering the multihadron cross section (fTOt . and the accuracy of 
the theoretic formulae employed for the calculation of the cross 
section. Using various initial values for the total and electron 
width we have verified that the iteration procedure converges 
fast introducing a negligible systematic uncertainty. To reduce 
the statistical error on the 4^(28) mass, the interference param- 
eter A was fixed in the fit at the value of 0.13 corresponding to 
Eq. (l22l) . Releasing the A parameter in the fit shifts the Ff^ x Si, 
value by -0.23%. This quantity can be used as an estimate of 
the influence of quasi-two-body decays with correlated inter- 
ference phases mentioned in Section 15.21 The accuracy of the 
resonance term in Eq. ( fTOl i is about 0.1% (Section |5]), and an- 
other 0.1% should be added because of the accuracy of radiative 
correction calculations in Ref. ifTsIl . The quadratic sum of these 
three contributions is about 0.3%. 

8. Inaccuracy due to interference in hadronic cross section 

The fits done with a floating interference parameter gave 

^ = 0.21 ±0.07 + 0.05 

The systematic uncertainty is mainly due to the beam energy 
determination and stability of the cross section measurement. 
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Table 4: Systematic uncertainties of the detection efficiency due to ip(2S ) decay simulation in % for three scans. The correlated parts 
of the uncertainties are also presented. The uncertainties for the first and second scans are assumed to be fully correlated. 



Source 


Scan 1 


Scan 2 


Scan 3 


Common 


Measured branchings 


0.4 


0.4 


0.3 


0.3 


JETSET ambiguities 


0.8 


0.8 


0.8 


0.8 


Track reconstruction 


0.4 


0.4 


0.7 


0.4 


Selection criteria 


0.3 


0.3 


0.3 


0.3 


MC statistics 


0.1 


0.1 


0.1 




Sum in quadrature 


s;1.0 


a=1.0 




a=1.0 



Table 5: The dominating systematic uncertainties in the Y^e x S/, product for three scans (%). The correlated parts of the uncertainties 
are also presented. The inaccuracy of about 0.9% due to possible interference phase correlation is not included. 



Source 


Scan 1 


Scan 2 


Scan 3 


Common 12 


Common 123 


Absolute luminosity measurements 


1.6 


1.7 


1.2 


1.6 


0.5 


ifr{2S) decay simulation 


1.0 


1.0 


1.1 


1.0 


1.0 


Detector response 


Trigger efficiency 


0.2 


0.2 


0.2 


0.2 


0.2 


Nuclear interaction 


0.2 


0.2 


0.3 


0.2 


0.2 


Cross talks in VD 


0.1 


0.17 


0.1 


0.1 


0.1 


Variation of cuts 


0.5 


0.3 


0.6 


0.3 


0.3 


Accelerator related effects 


Beam energy determination 


0.15 


0.18 


0.6 


0.15 


0.15 


Non-Gaussian energy distribution 


0.2 


0.2 


0.2 


0.2 


0.2 


Residual background 


<0.1 


<0.1 


<0.1 


<0.1 


<0.1 


Other uncertainties 


0.3 


0.3 


0.3 


0.3 


0.3 



Sum in quadrature ^2.0 ^2. 1 ^1.9 ^2.0 ^1.3 



This result does not contradict to the assumption of the uncor- 
rected interference phases (Section 15.21 ) but can still indicate 
the presence of some phase correlations. 

In order to evaluate a possible deviation of the f /, value from 
the sum of hadronic partial widths F/,, we performed a Monte 
Carlo simulation according to Eq. ( fTSl ) and ( fT9l l. The set of de- 
cay modes and the corresponding branching fractions were ob- 
tained using the event generator described in Section ISTl with 
the number of different decay modes exceeding one thousand. 
Damping of sin 0„, and cos (f)„ because of the averaging ( was 
ignored. 

The inaccuracy in F/, was estimated in the Bayesian ap- 
proach. Possible sets of interference phases 0„, are character- 
ized by the central value <p and the band width A<p. It was as- 
sumed that probabilities of all (p and A<p values are equal. The 
result is not sensitive to assumptions on the band shape. 

The Monte Carlo procedure consists of two steps. At the 
first step in each Monte Carlo sampling a set of cp,,, was gener- 
ated for random <p and A<p and the corresponding A value was 
calculated. The values of 0, and A were saved for the second 
step. The first-step distribution in A is symmetric relative to the 
most probable value of 0.13 corresponding to Eq. (fT4l i. At the 
second step an acception-rejection (Von Neumann's) method 



was employed to reproduce the A distribution matching the re- 
sults of the measurement: the Gaussian with the average value 
of 0.21 and the width of 0.086. The second-step distribution 
in (f;, - F/,)/r/, is peaked at zero. It is not Gaussian but 68% 
percent of the accepted samplings are contained in the interval 
of +0.009, thus we concluded that the inaccuracy due to pos- 
sible interference phase correlations in r^e x S/, value is about 
0.9%. Since the Bayesian approach was employed, appearance 
of the new information on the interference in the hadronic cross 
section can change this estimate. 

It should be noted that the inaccuracy estimated in this sec- 
tion is not specific for our results on the ^^^(25 ) partial and total 
widths but is shared by many results obtained in other exper- 
iments using the tfr(2S) multihadron cross section. The most 
precise of them is the result of CLEO on the i//{2S ) — > J/iffn^n:' 
branching fraction |37], its quoted accuracy is 2.2%. That con- 
cerns, in particular, the works (0,[3tl- 

9. Averaging of scan results 

The systematic errors on the F^c x S;, values for three scans 
and the estimates of their correlations are presented in Table |5] 
The correlation of errors is a difficult issue. In non-obvious 
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cases the most conservative approach was used assuming that 
the correlated part corresponds to the minimal uncertainty in 
scans for a given uncertainty source. 

To obtain the resulting value of the x Si, product, the 
scans were treated as independent experiments. The individual 
r^e X Si, values were weighted using their statistical errors and 
uncorrected parts of systematic errors. Such procedure takes 
into account the random behaviour of uncorrected systematic 
errors thus converting them to statistical. Correspondingly, the 
systematic errors of individual scans reduce to their common 
part. The formal weighting recipe for the parameter x Si, is 
given below: 

<r,,xs,,> = ^w,-(r,,xs,,)i , 

2 'V 2 2 

= 1/(0-L,,+C^syst,,-C^syst.o)' (32) 

where crsyst,o denotes a common part of systematic uncertain- 
ties. The recipe preserves the total error of the result. 



The statistical uncertainty is significantly reduced compared to 
that reached at VEPP-4M in 2002, the systematic one is approx- 
imately the same. Since the systematic and model errors are 
correlated, the combined analysis of the 2002, 2004 and 2006 
data has to be performed. It will be described in a dedicated 
paper the result of which should supersede the results presented 
above and in Ref. [10.1 . The reduction of the model dependence 
is expected. 
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10. Summary 



References 



The parameters of i//(2S) have been measured using the data 
collected with the KEDR detector at the VEPP-4M e+e" col- 
lider in 2004 and 2006. Our final result for the x Si, product 
is 



r,, X Si, = 2.233 + 0.015 + 0.037 ± 0.020 keV. 

The third error quoted is an estimate of the model dependence 
of the result due to assumptions on the interference effects in the 
cross section of the single-photon e^e" annihilation to hadrons 
explicitly considered in this work. Implicitly, the same assump- 
tions were employed to obtain the charmonium leptonic width 
and the absolute branching fractions in many experiments. This 
quantity was measured in several experiments but only the re- 
sult of MARK-I 113 811 . an order of magnitude less precise than 
ours, was pubUshed in such a form. Usually the Tg^ x Si, prod- 
uct is converted to the electron width value using existing re- 
sults on the branching fraction to hadrons Si, or the leptonic 
branching fractions. 

Using the world average values of the electron and hadron 
branching fractions from PDG ll24ll we obtained the electron 
partial width and the total width of i/j{2S): 



Tee = 2.282 ± 0.015 + 0.038 + 0.021 keV, 
r = 296 + 2 + 8 + 3 keV . 

These results are consistent with and more than two times more 
precise than any of the previous experiments. 

The result on the t//(2S) mass obtained in this work 



M : 



3686.1 14 ±0.007 ±0.011 to^l MeV. 
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